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ABSTRACT 

We present the first X-ray study of NGC6791, one of the oldest open clusters known (8 Gyr). Our Chandra 
observation is aimed at uncovering the population of close interacting binaries down to Lx ~ 1 x 10 3() erg s _1 
(0.3-7 keV). We detect 86 sources within 8' of the cluster center, including 59 inside the half-mass radius. We 
identify twenty sources with proper-motion cluster members, which are a mix of cataclysmic variables (CVs), 
active binaries (ABs), and binaries containing sub-subgiants. With follow-up optical spectroscopy we confirm 
the nature of one CV. We discover one new, X-ray variable candidate CV with Balmer and He II emission lines 
in its optical spectrum; this is the first X-ray-selected CV confirmed in an open cluster. The number of CVs 
per unit mass is consistent with the field, suggesting that the 3^1 CVs observed in NGC6791 are primordial. 
We compare the X-ray properties of NGC 6791 with those of a few old open (NGC6819, M67) and globular 
clusters (47 Tuc, NGC 6397). It is puzzling that the number of ABs brighter than 1 x 10 30 erg s _1 normalized 
by cluster mass is lower in NGC 6791 than in M 67 by a factor ~3-7. CVs, ABs, and sub-subgiants brighter 
than 1 x 10 3l) erg s are under-represented per unit mass in the globular clusters compared to the oldest open 
clusters, and this accounts for the lower total X-ray luminosity per unit mass of the former. This indicates that 
the net effect of dynamical encounters may be the destruction of even some of the hardest (i.e. X-ray-emitting) 
binaries. 

Subject headings: open clusters and associations: individual (NGC 6791); X-rays: binaries; stars: activity; 
binaries: close; cataclysmic variables 



1. INTRODUCTION 

X-ray emission of late-type stars arises in hot gas contained 
by magnetic loops above the stellar surface. These loops are 
produced by the interaction of convective motion and differ- 
ential rotation in the stellar envelope; hence the magnetic ac- 
tivity of a late-type main-sequence star — and with it the X- 
ray emission — is fo und to be highe r in rapidly rotating stars 
(see, for example, Pallavicini 1989). As stars age their rota- 
tion slows down, and it was a surprise that a ROSAT study of 
the old open cluster M 67 d iscovered a large number of X-ray 
sources (Bellon Tet al.|1993) . Comparison with optical obser- 
vations provided the explanation: in an old cluster stars may 
rotate rapidly when tidal forces sp in them up tow ards coro- 
tation with the binary revolution (Bello nTet al.|1998) . X-ray 
studies of old open clusters help in identifying such, otherwise 
inconspicuous, active binaries (ABs) in which tidal forces are 
or have been effective, and contribute to the understanding of 
tidal interaction. Since the typical timescale for tidal interac- 
tion rapidly incr eases w ith the ratio of the orbital separation 
to stellar radius (Zahn 1989), only stars in binaries with rel- 
atively short periods or relatively large radii experience the 
effects of tidal locking. X-ray observations of ABs in old 
clusters thus probe the populations of hard binaries, and the 
relative number of such systems found in clusters with dif- 
ferent properties is an important tool for understanding the 
effects of the cluster environment on the binary content. 

Thanks to prol onged optical studies, e.g. in the WIYN Open 
Cluster Study (Mathieu 2000), our knowledge of the binary 
population in clusters is reaching completeness levels which 



enable a detailed comparison with theoretical models of the 
evolution of stellar clusters and of the binary population in 
them. To explain the large number of blue stragglers ob- 
served in M 67 (4 Gyr), the model by |Hurley et aLl ((2005) 
requires a large initial number of short-period binaries that 
act as seeds from which blue stragglers are formed. As the 
simulated hard-binary fraction and period distribution do not 
evolve much, the close bi naries rem ain abundant throughout 
the life time of the cluster. Geller & Mathieu ( 2012 ) point out 
that comparison with the observed binaries in the old cluster 
NGC 188 (6.5 Gyr; |Meibom et al.|20 09) finds no evidence of 
such a high fraction of close binaries. The number of clus- 
ters studied is still small, as is the number of detailed models, 
and more study is required be fore a final ve rdict on these and 
other discrepancies (see Geller & Mathieu ( 2012| > for details) 
can be made. 

By now, observations of globular clusters with the Chandra 
X-ray Observatory have uncovered hundreds of faint (Lx < 
10 33 erg s _1 ) close binaries, only part of which have been clas- 
sified so far. As these binaries are an important driver of the 
long-term evolution of a cluster, understanding the properties 
and frequency of the various binary classes is a major goal 
of investigating (globular) cluster X-ray sources. Sometimes 
the outcome of a dynamical encounter is a binary or multi- 
ple system in a very unusual configuration, and this allows 
us to study source types that are not, or only rarely, found in 
the field. Besides ABs, which are the dominant X-ray source 
class in old open clusters, globular-cluster X-ray sources in- 
clude cataclysmic variables (CVs), low-mass X-ray binaries 
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in quiescence (qLMXBs), and milli-second pulsars (MSPs). 
The correlation of the numbers of qLMXBs and CVs with the 
stellar encounter frequency reveals the important role th at dy- 
namical encounters play in the creati on of such binaries ( |Poo-| 
ley & Hut|2006l |Heinke et al.|2006) , but the net balanceBe^ 



tween formation and destruction is less clear. For a large part 
this is due to the lack of a reference point that quantifies their 
number densities in a dynamically inactive environment. Ide- 
ally, the comparison should be made against the Galactic field 
where stellar densities are low, but this is complicated because 
of the generally limited information on distance and age for 
stars in the field, and the intrinsic scarcity of sources such 
as qLMXBs. Old open clusters provide an alternative that is 
worthwhile to explore, but so far only few have been stud- 
ied in detail in X-rays. Moreover, it should be kept in mind 
that dynamical encounters cannot be completely ignored in 
open clusters. Observational indications for the occurrence 
of dynamical interactions have come from individual systems 
whose properties are difficult to explain by binary evolution as 



it would take pla ce outside a cluster environment (e.g. |van den 
Berg et al. 2001 ) 



In this paper we present the results of the first X-ray ob- 
servation of NGC6791, which, at an age of about 8 Gyr, is 
one of the oldest known open clusters in our Galaxy. Be- 
sides the Chandra data, we obtained optical spectra to classify 
candidate optical counterparts to the detected X-ray sources. 
NGC6791 has been studied extensively in the optical. The 
rich body of available literature has proven to be very useful 
for the classification of our Chandra sources, and for sepa- 
rating cluster members from non-members. For the study of 
interacting binaries it is a promising target: before this work 
it was known to harbor two of only thr ee spectroscopically - 
confirmed CVs found in open clusters (Kaluzny et al. 1997), 
and many optical-variability studies have reveale d dozens of 
close binarie s in the field of the NGC6791 (see de Marchi 
et al.| ( [2007] ) and references therein). It is therefore interest- 
ing to investigate if the X-ray source content of NGC6791 is 
as rich and varied as was found for the few old open clusters 



studied previously, v iz. M67 ( |Belloni et aTp9"93l [19981 [van 
den Berg et al.|2004| > and NGC 188 ( |Belloni et al.|1998||Gon 



d. 



doin|2005| l. The cluster lies at a distance of ~4. 0-4.3 kpc and 



is reddened by E(B - y)=0.09-0.16; see e.g. Carraro et al. 



(2006),|Basu et aLlpOTT), and Brogaard~eTaT1(|201 1|> forre^ 



cent determinations of the cluster parameters. Throughout the 
paper we adopt a distance of 4. 1 kpc and a constant reddening 
E(B - V) = 0.14, which corresponds to a neutral-hydrogen 
column density N H = 7.8 x 10*° cnT 2 (|Predehl & Schmitt| 



1995). The reddening could vary across th e area of the cluster 
(Platai s et al.|201 l[|Brogaard et al.|2012| > but the effect is too 



small to have a significant impact on our results. NGC 6791 
stands out in several ways, including its high age and large 
mass. Remarkably for such an old cluster, the metallicity is 
higher than solar ([Fe/H]x!+0.4). 

The X-ray and optical observations and data reduction are 
described in Sect. [2] followed by the data analysis in Sect. [3] 
We detected various types of sources in the cluster which are 
presented in Sect. |4] Our results are discussed in the context 
of the populations of interacting binaries in star clusters in 
Sect. |5] We finish with our conclusions in Sect. [6] 



2. OBSERVATIONS AND DATA REDUCTION 
2.1. X-ray observations 




Figure 1. V-band image of NGC 679 1 from Stetson et al. 1 2003 1 showing in 
black the outline of the ACIS chips (labelled) as in our Chandra o bservation. 
The small wh ite, innermost circle indicates the cluster center from |de Marchi| 
|et al.| j2007) . The larger white circles mark the area inside the half-mass 
radius r/, (4:42), and the area inside 8', which was analyzed in this paper. 
North is up, east to the left. 

We observed NGC 6791 with the Advanced CCD Imaging 
Spectrometer (ACISflon Chandra from 2004 July 1 20:51 
UTC until July 2 10:49 UTC for a total exposure time of 48.2 
ks (ObsID 4510). The observation was done in very faint, 
timed exposure mode with a single-frame exposure time of 
3.2 s. In order to achieve optimal sensitivity below ~1 keV, 
we placed the central part of the cluster on the backside- 
illuminated S3 chip; the neighboring SI, S2, S4, 12, and 
13 chips were used to cover the outer parts of the clus- 
ter. The cluster center l ies at approximately a=19 h 20 m 53 s , 
(5=+37°46'18" (J2000; [deMarchi et al.|2007] l; we shifted the 
center Q'.5 away from the S3 aimpoint in the —Y (approxi- 
mately south-east) direction so that a larger part of the cluster 
could be imaged on a single chip. Platais et al. (201 1 ) derive 
a cluster half-mass radius r/, of 4'42 + i)'.{)2. Our observation 
covers almost the entire area inside r/, (> 98.5%) with either 
the S 3 or the S 2 chip; only the southernmost tip falls off the 
detector area. Fig.fTlshows the outline of the ACIS chips over- 
layed on an opticalimage of the cluster. 

Our data reduction started with the level- 1 event file pro- 
duced by version 7.6.7.2 of the Chandra X-ray Center pro- 
cessing pipeline. We used the CIAO 3.4 package with the 
CALDB 3.3.0 calibration files for further processing follow- 
ing standard procedure^] including the very-faint mode back- 
ground cleaning. In order to improve the precision of the 
source positions we removed the randomization of event co- 
ordinates that is applied in standard processing. Towards the 
end of the observation, the background level between 0.3 and 
7.0 keV increased by a factor of ~5 as a result of a background 
flare. After removing this interval from the observation, the 
exposure time is reduced to 42.9 ks. 

2.2. Optical spectroscopy 

We obtained low-resolution spectra of candidate optical 
counterparts to guide the classification of the X-ray sources. 

1 http://cxc.harvard.edu/proposer/POG/html/ACIS.html 

2 http://cxc.cfa.harvard.edU/ciao3.4/guides/acis_data.html 
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A total of 16 candidate counterparts brighter than V « 18.3 
were observed with the FAST long-slit spectrograph on the 
1 .5-m Tillinghast telescope on Mt. Hopkins on 9 nights be- 
tween 2005 June 7 to September 2. We used the 300 lines 
mm -1 grating, resulting in a wavelength coverage from 3480 
to 7400 A and a 3-A resolution. Exposure time was chosen 
to achieve a signal-to-noise ratio S /N > 20 for V < 17.5, and 
10 < S /N < 20 for fainter sources. FAST spectra were ex- 
tracted and wavelength-calibrated with a dedicated reduction 
pipelin^] Flux standards were observed on the same nights 
as the science targets. 

Candidate optical counterparts fainter than V ~ 17 were 
observed with the fiber-fed multi-object spectrograph Hec- 
tospec on the 6.5-m Multi-Mirror Telescope. Use of the 270 
lines mirT 1 grating resulted in spectra that cover 3700 to 9150 
A with a resolution of 6 A. A total of 16 candidate counter- 
parts were observed on the nights of 2005 May 13 and July 
4-6. Each setup was repeated 4 to 5 times with individual 
exposures of 900 s. Exposures were combined, and the spec- 
tra extracted and wavelength-calibrated, with the Hectospec 
reduction pipelin^] Sky background spectra for each setup 
were created by combining the spectra of fibers positioned 
at off-source locations. Hectospec observations did not in- 
clude flux standards. To correct for the instrumental response 
we constructed sensitivity curves using the spectra of the sub- 
dwarf B (sdB) st ars B 4, B 5 and B 3 in NGC 6791 ( |Kaluzny| 
& Udalski|1992| l, included in the May 13, July 4 and July 5- 
6 setups, respectively. We assumed that their intrinsic spec- 
tra can be described by blackbodie s with tempe r atures as de- 
termined for these sdB stars by Liebe rt et al.| ([T994). The 
blackbody spectra we re then no rmalized to reproduce the ob- 
served V magnitudes ( Stetson et al. 2003) after correction for 
the cluster reddening. This way we achieved an approximate 
absolute-flux calibration of the target spectra after applying 
the resulting sensitivity curves. 

Flux-calibrated target spectra were assigned spectral types 
by comparison to standard spectra with similar resolution (e.g. 
Jacoby et al.|1984| >. 



3. ANALYSIS 

3.1. X-ray source detection and extraction 

We restrict the analysis to the area within 8' of the cluster 
center. Sources further away have relatively large positional 
errors, which complicates the optical identification. We fo- 
cus on the area inside r/, which is close to the largest circular 
area around the cluster center that is entirely covered by the 
observation. 

We performed source detection in a broad (0.3-7.0 keV), 
soft (0.3-2.0 keV) and hard (2.0-7 keV) energy band, also 
used in our Chandra study of M 67 ( van den Berg et al. 2004 ), 
to facilitate comparison. The CIAO detection routine wavde- 
tect was run for scales of 1.0 to 1 1.3 pixels, in steps increasing 
by a factor V2, with the larger scales appropriate for large off- 
axis angles where the point-spread function (PSF) becomes 
significantly broader. We computed exposure maps for the 
response at 1 keV to account for spatial variations of the sen- 
sitivity. The wavdetect detection threshold was set to 10~ 6 , 
from which we expect two spurious detections per detection 
scale (so sixteen spurious detections in total) in the area that 
we consider here. Wavdetect positional errors reflect the sta- 



tistical uncertainty in centroiding the spatial distribution of 
the detected events of a given source, but do not include sys- 
tematic errors that result from data processing and PSF asym- 
metries at large offset angles. To c ompute positional uncer- 
tainties we therefore adopt Eq. 5 in |Hong et al.| (f2005). This 
formula relates the 95% confidence radius on the source posi- 
tion to the wavdetect counts and offset angle from the aim- 
point, and is based on extensive simulated detections of arti- 
ficial sources. Combination of the broad, soft, and hard-band 
source lists results in a master catalog of 86 distinct sources 
within 8' of the cluster center, of which 59 lie inside r/,. Ta- 
ble [TJ summarizes their basic properties. To investigate the 
validity of the sources, we also ran wavdetect with a thresh- 
old of 10~ 7 or an expected number of spurious sources of 1.6. 
The fourteen sources not detected in this run are marked with 
an asterisk in Table [T] In this paper we adopt a short version 
of the source names to refer to the sources instead of using 
their official CXO names; see columns 1 and 2 in Table [T| 
The shorthand names are assigned by first sorting the sources 
within r/, on net counts, and then the sources between r/ 2 and 



We used the acis ^extract package (version 3.107.2; |Broos| 
|et al.|[2002] > to measure net source counts. Events between 
0.3 and 7.0 keV were extracted from a region corresponding 
to 90% of the PSF at 1.5 keV; for a few sources this region 
was reduced to avoid contamination by a close neighbor. The 
background was determined from a source-free annulus cen- 
tered on the source position. We convert net count rates to 
fluxes within Sherpajwith arf and rmf response files appro- 
priate for the chip location and source-extraction area of each 
source. Column 7 of Table [T] lists absorption-corrected («) 
fluxes Fx,u computed under the assumption that the under- 
lying spectrum is an optically thin plasma (described by the 
xsmekal model) with kT = 2 keV. For source CX 29 that lies 
near the aimpoint 1 count s corresponds to 6.8 x 10~ 12 erg 
cirT 2 s . This temperature is ap propriate for the most ac - 



3 http://tdc-www.harvard.edu/instruments/fast 

4 http://tdc-www.harvard.edu/instruments/hectospec/ 



tive ABs in the cluster (see e.g. |van den Berg et aL 2004), 
while too high for the least active ABs and too low for 
most CVs. For a 1-keV Mekal model, a power-law spec- 
trum (powlawld) with photon index T — 2, and a 10-keV 
thermal-bremsstrahlung model (xsbremss) the conversion fac- 
tor is 24% smaller, and 30% and 42% larger, respectively. 
In each case we account for a column density equal to the 
cluster value using the xsphabs model. For the adopted dis- 
tance to NGC 6791 of 4.1 kpc, a 3-count detection limit cor- 
responds to a limiting unabsorbed X-ray luminosity Lx,u of 
(0.7 - 1.4) x 10 30 erg s -1 cirT 2 (0.3-7 keV) where the range 
corresponds to the choice of models specified above. 

Acis-extract performs a Kolmogorov-Smirnov test on the 
event arrival times to test for variability. Two sources are thus 
found to be potentially variable, with probabilities that their 
count rates are constant (p k-s) smaller than 2.5%: CX 19 
(a candidate CV, see Sect. 4.2| > for which 15 of a total of 
19 counts are detected within 5.5 hours (/?k-s = 0.018), 
and CX37 (no optical counterpart found) for which all 9 
counts are detected in the first 7 hours of the observation 
(Pk-s = 0.024). 

3.2. X-ray spectral properties 

Only CX 1 has sufficient counts to allow a constraining 
spectral fit. A spectrum was extracted from the acis-extract 

5 http://cxc.harvard.edu/sherpa 
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Table 1 

Chandra sources within 8' of the center of NGC 679 1 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


cx 


CXOUJ 


a(J2000) 


<5(J2000) 




Offset 


Counts 




£50 


Opt 






(deg) 


(deg) 


(") 


(') 




(erg citT 2 s~') 


(keV) 










Sources inside the half 


-mass 


radius ry, 








1 


192044.9+374640 


290.187175 


+37.777870 


0.32 


1.64 


213 ± 15 


318.0 


1.3 + 0.1 


+ 


2 


192039.8+374354 


290.165914 


+37.731752 


0.39 


3.54 


151 ± 13 


242.9 


1.00 + 0.03 


+ 


3 


192035.7+374452 


290.148996 


+37.747852 


0.39 


3.70 


125 ± 12 


260.1 


1.4 + 0.1 


+ 


4 


192056.3+374613 


290.234628 


+37.770331 


0.34 


0.66 


1 16 ± 11 


170.7 


1.4 + 0.1 


+ 


5 


192047.3+374318 


290.197161 


+37.721768 


0.43 


3.20 


72 ±9 


118.0 


1.4 + 0.1 




6 


192114.4+374530 


290.310342 


+37.758603 


0.63 


4.32 


72 ±9 


112.0 


1.4 + 0.1 




7 


192052.3+374550 


290.218110 


+37.764148 


0.35 


0.47 


53 ± 8 


77.9 


1.2 + 0.1 


+ 


8 


192038.2+374441 


290.159518 


+37.744838 


0.46 


3.33 


48±7 b 


221.7 


1.6 + 0.2 


+ 


9 


192058.4+375008 


290.243732 


+37.835590 


0.56 


3.98 


37 ± 7 


59.8 


1.6 + 0.3 


+ 


10 


192037.3+374612 


290.155726 


+37.770073 


0.43 


3.09 


34 ±6 


69.3 


1.4 + 0.2 





Note. — Columns: 1) source number; 2) source name; 3 and 4) source coordinates in decimal degrees after applying the boresight correction of Ao- = 
-(X'06 ± (X'06, = -0721 ± ff'04 (Chandra minus optical); 5) 95% uncertainty radius on the source position; 6) angular offset from the cluster center; 7) net 
counts in the 0.3-7 keV band; 8) unabsorbed flux (x 10~ 16 erg s cirT 2 ) in the 0.3-7 keV band for the assumption that the source spectrum is a 2-keV Mekal 
plasma seen through a neutral-hydrogen column of density_7.8 X 10 20 cm -2 ; 9) median energy in the 0.3-7 keV band (only for sources with >5 counts); 10) flag 
for the detection of a candidate optical counterpart. TableTTlis available in its entirety in a machine-readable form in the online journal. A portion is shown here 
for guidance regarding its form and content. 



source region with the CIAO tool psextract, and was grouped 
to have at least 20 counts bin -1 to warrant use of the^ 2 statis- 
tic; the background contribution (<1 count) can be ignored. 
Given the classification of CX 1 as an activ e gal actic nucleus 
(AGN) based on its optical spectrum (Sect. 4.6 1, we try to fit 
the data with an absorbed power-law and find an acceptable 
result for a photon index F = 1.9 + 0.3 and a column den- 
sity jV h < 1.2 x 10 21 cirT 2 (1-cr upper limit; x 1 = 9.5, 6 de- 
grees of freedom). The limit on Nu is consistent with the inte- 
grated Galactic column density in the direction of NGC 6791 
(9 x 10 20 cnr 2 , Sch legel et al.|l998|> whil e the value for T is 
typical for an AGN (e.g. Tozzi et a l.|2006) . 

The remaining sources have too few counts for useful con- 
straints by spectral fitting. Instead we use the method of quan- 
tile analysis developed by Hong et al. (2004), where the me- 
dian energy (£50) and the 25% and 75% energy quartiles (£25 
and £75) of the source photons are used to characterize spec- 
tral properties. The advantage of using energy quantiles as 
opposed to comparing counts in pre-defined energy bands by 
means of hardness ratios is that the errors on the diagnostics 
are less sensitive to the underlying spectral shape 

3.3. Optical identification 
3.3.1. Cross-identification against the Stetson catalog 

We looked for optical counterparts in the deep BVI photo- 
metric catalog of NGC 6791 compiled by Stetson et al. (2003; 
S03 hereafter), which covers the entire area studied here. The 
limiting magnitude is V ~ 24 for the central area but the 
sensitivity is lower for the outer regions. In an attempt to 
eliminate artifacts from the catalog, we removed entries with 
photometric -quality indicators that flag them as suspicious, 
viz. sources with |sharpness| > 1, and separation index < 0; 
see S03 for an explanation of these indicators and a motiva- 
tion for these criteria. 

We first measure and correct for the boresight, i.e. a possi- 
ble systematic offset between the astrometric reference frames 
of the Chandra and optical positions. The optical posi- 
tions are tied to the International Celestial Reference System 
(ICRS) with an rms accuracy of about 0727 (S03), but the as- 



trometric reference frame of any given Chandra observation 
as a whole is aligned to the ICRS with a typical 076 accu- 
racy (90% uncertainty^. We measure the boresight correction 
using matches with optical variables only. This minimizes 
the number of chance alignments as the light curves of many 
variables reveal them as close binaries and therefore as plau- 
sible X-ray emitters. A list of variables was compiled from 



the studies of |Kaluzny & Rucinski| ([T9 93 ), Rucinsk i et al. 
(1996), Mochej sETet al.| (|2002[ |2003| |2005>, |Bruntt et aL 
|2003| , [Hartman et al.| ( |2005) , and |de Marchi et al.| (|2007 ); 
we only selected variables that we could identify with ob- 
jects in the S03 catalog so that we could use the S03 posi- 
tions for the boresighting. The boresight correction is com- 
puted iteratively. We looked for matches to X-ray sources 
with more than 10 net counts and rgs < I". For each source, 
the match radius is set to be the quadratic sum of the error 
on its X-ray position, and the typical 1-cr error on the opti- 
cal positions (0727) scaled to a 95% error radius assuming a 
2D gaussian error distribution. The boresight is set to be the 
weighted (with r 9 2 ) mean of the X-ray-optical offsets of all 
matches found. Next, the X-ray positions are corrected for 
the boresight, and the matching is repeated, now including 
the statistical error on the mean boresight in the match radius. 
This procedure quickly converges to a boresight correction of 
Aa = -0706 ± 0706, A6 = -0721 ± 0704 (Chandra minus op- 
tical) based on ten matches with optical variables. Matching 
the boresight-corrected X-ray catalog to the clean S03 cata- 
log results in 5 1 optical matches (including 26 variables) for 
47 sources, out of the 86 X-ray sources in total. 

The probability to find a match in the S03 catalog by chance 
depends on the projected star density (which decreases with 
distance from the cluster center r) and on the match radius 
(which typically increases with r because rgs increases). To 
estimate the expected number of spurious matches, we divide 
the cluster in a central area defined by r < ry,, and an outer 
area of r/, < r < 8'. The clean S03 catalog gives a projected 
density of 0.030 optical sources arcsec~ 2 and 0.0092 sources 



1 http://cxc.harvard.edu/cal/ASPECT/celmon/ 



A Chandra X-ray study of NGC 6791 



5 



arcsec~ 2 for the central and outer area, respectively. If we 
multiply this by the total area covered by the match circles of 
the X-ray sources, we find that the expected number of ran- 
dom matches is 5.9 (center), and 3.8 (outer annulus). This 
amounts to 15% and 29% of the matches found in the respec- 
tive regions. On the other hand, a similar calculation shows 
that all matches with variable stars are likely to be real, with 
the average number of chance alignments < 0.1 both for the 
central and outer regions. 

In order not to overlook any matches, we repeated the 
matching using the entire S03 catalog, and inspected the re- 
gions around the Chandra sources in the optical fits image 
from S03 to discard matches with image artifacts or dubi- 
ous detections. We thus found nine extra candidate coun- 
terparts for eight Chandra sources, which are not included 
in the clean catalog because the values of the quality flags 
slightly exceeded the adopted cutoff limits, the object is re- 
ally extended, or because the object is faint and lies close to a 
relatively bright star. Since selection of these additional coun- 
terparts was not done in a very systematic way but relies on 
visual inspection, it is not trivial to estimate the number of 
random coincidences among the new matches. Two new can- 
didate counterparts are matched to a single Chandra source 
outside r/,, so at least one of the two must be spurious. The 
other seven are uniquely matched to seven Chandra sources 
inside r/,. Based on the higher source density in the entire S03 
catalog compared to the clean S03 catalog, one would expect 
to find ~1.5 extra spurious matches in the central area. Keep- 
ing in mind that some of these optical "sources" are not real, 
we estimate that at most one or two of the seven new matches 
inside ry, are spurious. 

In total the S03 optical survey provides 60 astrometric 
matches for 55 sources. 



3.3.2. HST imaging 
Bedin et al. (|2006)l used Hubble Space Telescope (HST) 



multi-epoch imaging with the Wide Field Channel on the Ad- 
vanced Camera for Surveys (ACS) to measure proper mo- 
tions of objects in a central 3'A x 3'A region of NGC 6791. 
We make use of their results to establish cluster member- 
ship o f ca ndidate counte rparts included in th e ACS images 
(Sect. |4~T) , and refer to |Bedin et aT] ( |2006| ) for details of 
the data and the proper-motion analysis. We also use these 
deep images to look for additional faint candidate counter- 
parts. The data set consists of F606W and F814W images 
taken on 2003 July 16 and 17 (GO-9815), and on 2005 Jul 13 
(GO-10471). An astrometric reference frame was created by 
combining the GO-9815 F814W images with the STScI Mul- 
tidrizzle software, which removes artifacts like cosmic rays 
and bad pixels, and corrects for the geometric distortion of 
ACS images. In order to align the coordinate system to the 
ICRS, we computed an astrometric solution based on the posi- 
tions of 273 unsaturated stars from the S03 catalog that could 
be identified in the stacked image. Fitting for zero point, plate 
scale, and rotation resulted in a solution where the rms residu- 
als of individual stars are 07035 in right ascension and 07032 
in declination; this is negligible compared to the X-ray posi- 
tional errors. The boresighted X-ray positions and error cir- 
cles were overlayed on this image, and photometry and proper 
motions for the astrometric matches were extracted from the 
source catalog created from the entire HST data set. We find 
seven new candidate counterparts for six X-ray sources. The 
new matches are all likely extra-galactic, given their proper 
motion or extended morphology. From the fact that two of 




Figure 2. ACS/WFC finding charts, each measuring 3'.'4 X 3'.'4 in size, made 
from the stacked GO-9815 F814W image representing 7024 s of exposure 
time in total. The combined X-ray/optical 95% error circles are shown as 
circles centered around the boresighted X-ray positions, while tick marks 
indicate the optical matches. Except for the bright candidate counterpart to 
CX 2 1 on the edge of the error circle, the se matches are found in the HST data 
only, and not in the Stetson et al. ( 2003 1 catalog. North is up, east to the left. 

the seven HST candidate counterparts are matched to Chandra 
sources with another candidate counterpart, we estimate that 
an appreciable fraction of HST matches could be random, on 
the order of 2/7 (^29%) or even more. As the astrometrically- 
calibrated image is not readily available, we show the find- 
ing charts of these additional identifications in Fig. [2] This 
brings the final tally of the optical identification to 53 candi- 
date counterparts for 48 of the 57 X-ray sources inside r/,, and 
14 candidate counterparts for 12 of the 29 sources between 
r/, and 8'. The properties of the candidate counterparts are 
summarized in Tables [2] and [3] and their locations in the opti- 
cal color-magnitude diagrams (CMDs) are shown in Figures|3] 
and [4] More information on the light-curve properties of the 
optical variables (column 10 of Tabled can be found by con- 
sulting the original references. Variables with names listed in 



t he for mat nnnnn_ * were first d iscovered by de Marchi et al. 
( |2007[ ); Table 1 in de Marchi et al. ( |2007| l gives an overview 
of the discovery papers for variables with names starting with 
a 'V or 'B', although updated light curves can sometimes be 
found in more recent papers mentioned in that table. 

3.3.3. Sources with multiple counterparts 

For seven X-ray sources we find two candidate counter- 
parts. In four cases— CX 27, CX57, CX78, CX84— one of 
the two is an optical variable, and given the low probability 
for a chance coincidence we consider the variables to be the 
true matches. The variables all lie closer to the X-ray source 
than the alternative match. We also tentatively list the clos- 
est match to the other three X-ray sources as the likely coun- 
terpart in Table [2] but more information is needed to firmly 
establish which, if any, of the two is the true match. For 
completeness we list the properties of the alternative identi- 
fications in Table [3] CX 34 is likely an extra-galactic source, 
as both HST counterparts are extended. The closest match to 
CX21 is a faint point source that is unrelated to the cluster; 
its HST proper motion is typical for a background galaxy. The 
other match is a proper-motion member on the main sequence 
(see also Fig. [2). CX82 matches with two bright (V < 15) 
stars. The nearest is a proper-motion non-member and a likely 
foreground K dwarf. 

3.4. Source classification 
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Figure 3. (V, B-V) and (V, V - I) color-magnitude diagrams of NGC 6791 with the candidate optical counterparts that are associated, or possibly associated, 
with the cluster marked in red. Candidate counterparts that are proper-motion members are indicated with large filled symbols, and sources without membership 
information with large open symbols or crosses: diamonds are (candidate) CVs, circles are (candidate) ABs, triangles are sub-subgiants, squares are stars without 
any indication of binarity, and crosses mark unclassified faint sources away from the main sequence. Stars in the field of NGC 6791 are plotted as dots, while 
stars with a proper-motion membership probability > 50 % (K. Cudworth, priv ate communication) are shown as small filled black circles. The insets zoom in on 
the crowded regions of the diagrams. Photometry is from Stetson et al. 1 2003 1. See the electronic edition of the Journal for a color version of this figure. 



Table 3 

Alternative optical counterparts 



CX 


OID 


dox 
(") 


V 


B-V 


V -I 


comment 


21 


11098 


0.78 


18.48 


0.93 


0.97 


Pk = 90 


27 


11365 


0.68 


16.58 


0.76 


0.94 


Pn = 0, F star 


34 




0.68 








extended 


57 


7395 


0.76 


21.02 


1.47 


1.96 


non-member, K star 


78 


569 


3.12 


16.50 


0.95 


1.03 


/V = 


82 


13395 


0.70 


14.66 


0.96 


0.98 




84 


15442 


3.01 


20.64 




1.7 





Note. — See Table[2]for the closest astrometric optical match to the X-ray 
source, and for the meaning of the columns. 



Many factors contribute to the classification of the X-ray 
sources. An extended morphology of the candidate opti- 
cal counterpart separates extra-galactic from galactic sources. 
Based on optical spectra one can immediately distinguish 
between AGN, accreting binaries with substantial accretion 
disks, and ordinary stars. Period and shape of the optical 



light curve constrain the nature and period of a close bi- 
nary. Proper-motion information, and to some extent opti- 
cal colors, can establish cluster membership and puts limits 
on the distance. To distinguish between cluster members and 
non-members we m ainly rely on the proper-motion study by 
Platais et al. (2011), which includes 41 of our candidate coun- 
terparts. For nine additional candidate counterparts member- 
ship could be esta blished from the HST study by Bedin et 
al. (2006, see Sect. 3.3.2 1 and from the results of the proper- 
motion study by K. Cudworth (private communication). 

We use the energy quantiles to constrain the underlying X- 
ray spectrum for sources with more than 15 net counts (0.3-7 
keV). For fainter sources the errors on the quantiles are too 
large to meaningfully say anything about their spectra. Quan- 
tile color-color diagrams (QCCDs) are shown in Fig. [5] Each 
panel shows two types of model grids, which indicate the ex- 
pected locations of a source in the QCCD for emission from 
a Mekal plasma (blue/yellow) and for a power-law spectrum 
(black/gray). The former is appropriate for the emission from 
hot coronae around active stars for which the plasma temper- 
ature kT can reach up to several keV. CVs can be anywhere 
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Figure 4. Color-magnitude diagram of NGC 679 1 with candidate counter- 
parts that are likely not associated with the cluster marked in red. Downward- 
pointing triangles are background galaxies, circles are (candidate) ABs, and 
squares are likely foreground stars without any indication of binarity. HST 
counterparts, for which photometry in V and / is not available, are not shown. 
See the electronic edition of the Journal for a color version of this figure. 

between soft or very hard; th e spectrum of a typ ical dwarf 
nova has kT w 2 - 10 keV (e.g. Byckling et al. 2010|l. Power- 



law spectra are appropriate for the harder emission from AGN 
(F = 1 -2). Grids are computed with Sherpa using the energy 
response of the source closest to the aimpoint. Variations in 
the shape of the grids as function of chip position are typically 
smaller than the errors on the measurements, and are ignored. 

The ratio of the fluxes in the X-ray and optical (V) band, 
or limits thereon, can help to classify a source, even in the 
absence of an optical match down to the limit of the obser- 
vations. We compute this ratio as follows: log(f x/F opt ) u = 
log(F x )u + (V - A v )/2 5 + 5.44. The zer opoint for the V- 
band flux is taken from Bessel ^t~aL"1 ( |1998] l. The total extinc- 
tion Ay was assumed to be equal to the cluster value, which 
underestimates (overestimates) the flux ratio for foreground 
(background) objects. Typically, stars and active binaries 
have \og(Fxl Fy) u < — 1, while AGN, CVs or other accreting 
binaries with unevolved late-type companions, and very ac- 
tive late-type dwarfs can have higher ratios (e.g. Stocke et al. 
[1991) . 

4. RESULTS 

4.1. Cluster members versus non-members 

Twenty candidate counterparts are proper motion members. 
Of these, all but one (CX 9) show signs of binarity. On this 
basis we consider at least 19 of them as the true counterparts 
to the Chandra sources. We derive the binary status and type 
of binary from the optical spectra, the optical light curves, or a 
location in the CMD on the binary main sequence, indicating 
that an unresolved multiple system is responsible for the de- 
tected light. The cluster binaries are a mix of CVs, ABs, and 
binaries below the sub -gian t branch, and are discussed further 
in Sects. 



4.2 4.3| and |4.4| respectively. CX9 is discussed in 
Cluster members are marked with filled symbols in 



Sect. 4.5 

theCMCs of Fig. [3] 



A total of 30 sources are not, or likely not, associated with 
the cluster. These are stars and binaries with proper motions 
that clearly set them apart from the members, background 
galaxies, and stars whose very red colors classify them as un- 
likely members . They are s how n in Fig. [4] and briefly dis- 
cussed in Sects. |4~5)|4 6| and |4.7| 

For ten candidate counterparts membership information is 
lacking or inconclusive. These sources include (candidate) 
ABs, a new candidate CV, and faint optical sources that lie 
to th e blu e of the main sequence; the latter are discussed in 
Sect. 4.7 Open symbols and crosses mark them in Fig. [3] 

In the following we discuss the sources by type of X-ray 
emitter. Unless stated otherwise, X-ray luminosities quoted 
in the text refer to the 0.3-7 keV band. 

4.2. Cataclysmic variables 

We have detected three CVs that belong to the cluster and 
discovered one CV candidate without membership informa- 
tion. For all four, the X-ray colors and luminosities are in the 
expected range for CVs. While the X-ray-to-optical flux ra- 
tios in Table 2] are also typical for CVs, their values can be 
misleading: the X-ray and optical data are not contemporane- 
ous, while CVs may show large variations in brightness on a 
time scale of weeks to months. 

CX4 is matched to the faint cluster member and optical 
variable 06289-9 that lies to the blue of the cluster main se- 
quence. Based on its op tical colors and t he detection of an 
outburst-like event, de Marchi et al. ( 2007[ ) already suggested 
that this star is a dwarf nova. To our knowledge the Hectospec 
spectrum in Fig. [6] which shows the Balmer lines clearly in 
emission, is the first spectroscopic confirmation of its CV na- 
ture. The X-ray quantiles suggest that the X-ray emission 
arises in a kT « 4 - 10 keV plasma (Fig. |3J). Thus L x in 
Table [2] which was computed for the assump tion of a 2-keV 
plasma, can be up to ~40% too low (see Sect. |3.2[ ). 

The new candidate CV is the X-ray variable CX 19, which 
is matched to a faint blue object with a Hectospec spectrum 
that shows Hell 4686 A and H/3 emission lines (Fig. [H]), 
and hints of Ha and H6 in emission. The X-ray luminosity 
(Lx,u ~ 6 x 10 30 erg s~'), the hard spectrum as suggested by 
its quantile values (Fig. BJ, and the high (Fx/Fy) ratio, are 
all consistent with a classification as CV. He II 4686 A is of 
comparable strength as H/3. This is seen in high mass-transfer 
rate systems (nova-like CVs), and in CVs containing white 
dwarfs with much stronger magnetic fields (i.e. >1 MG) than 
in dwarf novae systems. If we assume cluster membership 
we find an absolute magnitude M v m 10.2, which favors the 
explanation as a magnetic CV 

CX3 and CX 17 are th e known, spectros copically- 
confirmed CVs B8 and B7 (|Kaluzny et~aT| fT997]>. With 



Lx.u ~ 5.2 x I0 il erg s -1 , B8 is the brightest X-ray source 
in the cluster. The actual value of Lx is likely higher, as the 
X-ray quantiles suggest a plasma temperature (~4 keV) that 
is higher than our nominal value. B8 has been classified as 
a SUUMa dwarf nova based on the detection of several out- 
bursts of 1-2 mag in amplitude, an d a recent superou tburst 
(~3 mag) in the Kep ler light curve (Mochejska et al. 2003 
|Garnavich et al.||2011) . We include B8 among the cluster 
members. jPlatais et al.| ( [20TT| l estimated a 10% membership 
probability from proper motion, but closer inspection of the 
data reveals that a few detections of this star were contam- 
inated by a fainter close neighbor. When the affected data 
points are not considered, the proper motion of B8 is consis- 




log 10 (E 50 /0.3 keV) / log, (7 keV/0.3 keV) 



log, (E 50 /0.3 keV) / log, (7 keV/0.3 keV) 



Figure 5. Quantile color-color diagrams with model grids representing a Mekal plasma (blue/yellow) and a power-law spectrum (black/grey). For a certain 
choice of model, the plasma temperature or spectral index, and the column density can be estimated from the location of a source inside the grid. The median 
energy £50 can be read off from the top x-axis. Symbols are as in Figures [3] and [4] while black filled circles mark sources without optical counterparts. Due to 
their different energy responses, sources on the background-illuminated (Bl) S 3 chip and foreground-illuminated (FI) S 2 and S 4 chips are shown in separate 
panels. We include sources with more than 15 net counts (0.3-7 keV); error bars are shown only for the sources with the highest and lowest number of counts. 



tent with cluster membership. 

CX 17 or B7 is a solid cluster member. The X-ray quan- 
tiles indicate that kT » 2 keV. Most reports in the literature 
find B7 to be relatively bright with V ~ 18, which suggests 
a high mass-transfer rate; occasional small (0.5-1 mag) out- 
bursts, an d a drop in brightness o f ~3 mag in V have also 
been seen. Mochejsk a et al.| ( [2003| l propose that B7 is a nova- 
like variable of the VY Scl type, or perhaps a Z Cam dwarf 
nova mostly seen during periods of standstill between mini- 
mum and maximum brightness. Our Hectospec spectrum is 
qualitatively similar to the spectrum in |Kaluzny et al.| ( fl997| l, 
and likely taken when the mass-transfer rate was high and the 
disk optically thick: the spectrum shows broad Balmer ab- 
sorption lines with narrow emission cores, with the emission 
component dominant in Ho- (Fig. [6]). 

The X-ray-to-optical flux ratios and blue optical colors of 
CX18, CX25, CX28, CX36, CX65, and CX72 resemble 
those of CVs, and these sources are prime targets for optical 
follow-up spectroscopy. 

4.3. Active binaries 

We classify a total of twenty sources as likely or candidate 
ABs. Thanks to the extensive photometric variability stud- 
ies of NGC 6791, the orbital periods of many of these can be 
inferred from the light curves, either through eclipses, ellip- 
soidal variations, or spot activity. Among the cluster members 
we find eleven ABs, including eight with photometric periods 
(Table |2j. Most have colors that place them near the main 
sequence or, in the case of the W UMa contact binary CX 39, 
right at the turnoff. For CX 86, the 1-sigma errors on B—V and 
V - I are such that it could lie anywhere from just below the 
turnoff to the base of the red-giant branch. CX 23 lies securely 
above the sub-giant branch, and we discuss it separately be- 
low. NGC 6791 is at least 8 Gyr old; by comparison to the or- 
bital period versus eccentricity distribution of solar-type spec- 



troscopic binaries in the 6.5-Gyr old open cluster NGC 188 it 
is expected that, for main-sequence binaries, orbits up to at 
least 15 days have be en circularized at the age of NGC 6791 
(Mathieu et al. 2004). Therefore, since tidal synchroniz ation 
operates on a snort er time scale than tidal circularization ( |Hut| 
1 1981 1 jZahn 1989 ), main-sequence binaries with periods be- 
low 15 days are also expected to be tidally locked. All of our 
main-sequence ABs with measured phot ometric periods have 
periods shorter than 4.9 d; in Sect. 5.2 we discuss how this 



could be due to the period versus X-ray luminosity relation of 
ABs and our detection sensitivity. CX50, CX54, and CX79 
are not identified with known variables, but they lie on the bi- 
nary main sequence or have Ha in emission in their spectra, 
which is a signature of enhanced magnetic activity. We note 
that in the V versus B — V CMD, CX 54 lies slightly to the 
blue of the main sequence. This is not expected for the com- 
bined light of two main-sequence stars, so its classification as 
AB is less secure. Perhaps this star is a non-member after all, 
with a proper motion that is similar to that of the cluster stars. 
Radial-velocity measurements can provide more clarity. 

Stars that have evolved past the turnoff, i.e. have evolved 
to larger radii, can in principle circularize wider orbits than 
main-sequence stars within a given time span as the circu- 
larization time scale is very sensitive to radius. For the sub- 
giant AB CX23 the eccentricity is unknown, and estimates 
of the binary period are on ly known from photometric vari- 
ability BrunttetaL ( 2003] ) quote a period of ~ 12.5 d, while 
Mochejska et al. (2005) find nearly twice that value, ~23.9 d. 



From the V and B—V measurements of the optical counterpart, 
we estimate that the radius of the primary is at most ~2.8 R Q 
( |Flower|1996| ); as both components of the binary likely con- 
tribute to the measured optical light, the actual radius of the 
primary is smaller. If we now look at the diagnostic diagram 
in Verbunt & Phinney (1995, Fig. 5a; appropriate for a pri- 
mary mass of 1.25 M Q , i.e. somewhat larger than the turnoff 
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Figure 6. Hectospec spectra of t wo CVs and one c andidate CV. Top: our optical spectrum of CX 4 = 06289-9 confirms this optical variable and cluster member is 
a CV, as was already suggested by de Marchi et al. 1 2007). Middle: the spectrum of the known member CV CX 17 = B7 shows narrow emission lines superposed 
on a continuum with broad Balmer absorption features, indicating a high accretion rate. Bottom: the new candidate CV CX 19. To reduce noise, its spectrum has 
been smoothed with 3 pixels (A/1 a 3.6 A). Regions where imperfect sky subtraction left large residuals (around 5577A and 6300A) have been removed. 



mass in NGC 6791), it is clear that, to first approximation, 
the primary in CX23 is too small to have circularized a 23.9- 
d orbit. This means that in order to explain the X-rays, the 
rotation of the primary would have to be synchronized to the 
orbital motion around periastron, assuming the binary was not 
circular to begin with. However, this leaves the photometric 
period, which has been tentatively attributed to spot activity 
and thus reflects the primar y's rotation period, unexplained. 



m p i 

The shorter period found by|Bruntt et al.|(2~003|) is more con 



sistent with the detection of X-ray activity, as a 12.5 d orbit is 
short enough for tidal synchronization and circularization to 
have been achieved. The level of X-ray activity is still unusu- 
ally high for such a period (Sect. |5.2) . With radial-velocity 
measurements the exact period and eccentricity of the orbit 
can be verified. 

Three candidate counterparts without or with ambiguous 
membership information are classified as (candidate) ABs 
near the main sequence. CX56 is a 1.6-d spotted variable. 
If CX 49 is a member, it is almost certainly an AB owing to 
its location on the binary main sequence and Ha emission. 
CX58, for which we do not have an optical spectrum, lies 
on the cluster main sequence but could still be an AB in the 
cluster if the companion is faint. 

Six proper-motion non-members are classified as ABs. 
These include the W UMa contact binaries CX 24, CX 27, and 
CX78. The distance modulus inferred from the relation be- 
tween period, color, and absolute magnitude of contact bina- 
ries also place them in front of the cluster (Mochejska et al. 
[2003] l. 

4.4. Sub-subgiants 



Five sources are matched to proper-motion cluster mem- 
bers (p m - 99%) that stand out due to their optical photom- 
etry, which places them below or to the red of the sub-giant 
branch: CX 15, CX30, CX41, CX68, and CX77. Stars with 
similar properties have been labeled "sub-subgiants" or "red 
stragglers" owing to their unusual location in the CMD, which 
cannot be explained by the combined light of two ordinary 
cluster members. About two dozen sub-subgiants are known 
in open and globular clusters. Various mechanisms (past or 
ongoing mass transfer, dynamical encounters) have been sug- 
gested to explain their photometry but most systems still defy 
explanation. 

All five X-ray detected sub-subgiants in NGC 6791 are op- 
tical variables with periods between 3.2 and 13.8 d, and their 
puzzl i ng photometry has been pointed out before (Kaluzny 
|2003||Platais et al.|20TT) . CX 15 is matched to the eclipsing 
binary V9. Its striking light curve shows smooth, asymmetric 
brightness modulations that migrate in phase on a time scale 
of months; these have been attributed to large star spots, and 



indicate a h igh level of activity (see for example Mochejska 
|et al.|[2003] ). The other four systems show a single maxi- 
mum and minimum in the light curve when folded on the pe- 
riods cited in Table |2j and have been classified as plausible 
rotational spotted variable s (Mochejska et al.|2002||Kaluzny| 
|2003|[de Marchi et al.|2007) ; unambiguous proof of their bi- 
narity and orbital periods requires radial-velocity follow-up. 
While the photometry of CX41 can be accounted for by the 
light of three stars — in a bound multiple system or aligned by 
chance — the red colors of the other four require at least one 
anomalous star like an underluminous (sub)giant or an overlu- 
minous main sequence star. Their detection as X-ray sources 
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indicates ongoing interaction between companions in a close 
binary. The weak Ha emission in the optical spectra, and the 
relatively soft X-ray colors are more in line with coronal ac- 
tivity than with an accretion origin of the X-rays. Two more 
sub-subgiants in NGC6791 were pointed out by Platais et al. 
( I20TTT ). Star 13753 from S03 falls outside the Chandra field. 
Star 83 is not detected but it lies almost 9' from the observa- 
tion aimpoint so was observed with low sensitivity. 

4.5. Stars without indications ofbinarity 

CX9 matches with an early-K giant and cluster member. 
Our optical spectrum shows no indications of activity, nor is 
this star known to be an optical variable. Its X-ray colors are 
unusually hard for a coronal source (Fig.[3J). This giant is an 
interesting target for follow-up studies to investigate if it is 
a symbiotic-like interacting binary, in which the X-rays are 
powered by accretion. Alternatively, a faint quasar could be 
hidden under the point-spread function of this bright star, so a 
possible explanation is that the giant is a spurious match. 

Proper motions exclude cluster membership for the can- 
didate counterparts to CX7, CX 12, CX29, CX35, CX38, 
CX46, CX51, CX80, CX82, and CX84. These sources are 
plausible late-type field stars. The proposed counterparts for 
CX7, CX 12, and CX84 are reported to be long-term or ir- 
regular variables, but the light curves do not point at binarity; 
the remaining stars are not known to be variables. Our FAST 
spectra indicate that CX 35 and CX 46 are (sub)giants. Our as- 
sessment of luminosity class is mainly base d on the absence 
or pre sence of the CN 4216 A bandhead ( |Gray & Corbally 
2009 ). The log(fx/^v)u values, and the soft X-ray colors of 
CX7 and CX 12 agree with a coronal origin of the X-rays. 

The colors of the faint source CX 43 place it about 0.8 to the 
red of the lower main sequence and make it an unlikely cluster 
member. On the other hand, its proper motion is consistent 
with that of cluster stars. A radial-velocity based membership 
constraint can provide a more definite answer. Based on the 
red colors and high \og(Fx/Fy) u ratio, we suggest this source 
is a vey active late-type star. 

4.6. Background galaxies 

We classify twelve sources as extra-galactic on various 
grounds. The optical spectra of CX 1, CX 8, and CX 31 have 
broad emission lines typical for AGN, and we estimate the 
redshift z (see Table [2| b y comparison with th e quasar com- 
posite spectrum from Vanden Berk et al. ( 200 1 1 > . The counter- 
part to CX 45 has a redshifted stellar absorption-line spectrum 
and clearly is extended, even in the ground-based image from 
S03. The proper motion of the faint unresolved counterpart of 
CX21 is consistent with those of background galaxies in the 
field ( Bedin et al.|2006 1, but we note this source has a bright 
alternative counterpart on the edge of the error circle that is a 
cluster member (Sect. \333\ . CX 13, CX26, CX34, CX42, 
CX47, CX52, and CX53 look extended in the ACS images. 

Apart from CX 13 and CX34, which are too faint to get 
reliable optical photometry for, the extra-galactic sources 
lie away from the cluster main sequence and most have 
log{Fx/Fv) u Z -0.6. We note that the value of Nh used to 
calculate this ratio can be underestimated (and the flux ratio 
therefore overestimated) since extra-galactic sources can have 
significant intrinsic absorption. 

4.7. Uncertain classifications and unidentified sources 

Six sources are matched to faint, blue objects without mem- 
bership information: CX18, CX25, CX28, CX36, CX65, 



and CX72. As mentioned in Sect. 4.2 their optical colors 
and X-ray-to-optical flux ratios are compatible with those of 
CVs or other accreting binaries with low-mass companions, 
but also with AGN. Follow-up spectroscopy is needed to fur- 
ther investigate their classification. 

A total of 27 sources remain without candidate optical 
counterparts. For these sources we estimate a lower limit on 
the X-ray-to-optical flux ratio by assuming V > 24 (i.e. the 
limit of the S03 catalog), which gives log(Fx/Fv) u Z -0.9. 
This excludes a stellar coronal origin unless they are very ac- 
tive (flaring) K or M dwarfs, but leaves open the possibility 
of accreting compact binaries. An extra-galactic nature is the 
most likely option based on the expected number of AGNs. 
We compute this number using the cumulative sour ce den - 
sity versus flux (log N - log S ) curves from |Kim et al. ( 2007 1, 
which are derived from Chandra pointings at high Galactic 
latitudes. We assume a 3-count detection limit for the region 
within 3' from the aimpoint, a 5-count limit for larger off- 
sets, and an intrinsic power-law spectrum with photon index 
r = 1.7 that is absorbed by the integrated Galactic column 
density in the direction of the cluster (9 x 10 20 cm -2 ). This 
predicts ~56 AGN. For comparison, we detect 12 galaxies. 
This suggests that most of the unclassified and unidentified 
sources (33 in total) are AGN. 

5. DISCUSSION 

5.1. Cataclysmic variables in NGC6791 

NGC6791 is the open cluster where most CVs have been 
found so far. This may not be surprising because, of the open 
clusters that are old enough to find close binaries in a phase 
where they can be observed as CVs, NGC6791 is the most 
massive one. We check if primordial binaries can account 
for all observed systems, which comprise the three members 
(the dwarf novae B8 and CX4, and the nova-like variable or 
dwarf nova B7) and the new (possibly magnetic) CV candi- 
date CX 19. Recent observational estimates of the CV number 
density in the local field differ by a factor of a few but a re con- 
sistent within the (appreciable) errors: 4*5 x 10" 6 pc~ 3 



(|Preto-| 



rius & Knigge|2012 non-magnetic CV s only), 0.9^ X 10~ 5 

pc~ 3 ( [Rogel et al.|2008) , ~ 1 x 10~ 5 pc -3 ( |Grindlay et al.|2005| ). 
Adopting the local stellar mass density of 0.045 M Q pc from 
Robin et al. ( 2003| l, and scaling the resulting number of field 
CVs per u nit mass to the cluster mass (~5000-7000 M Q , see 
Sect. [5.4) , we expect 0.2-3.7 CVs in the cluster. Therefore 
the observed number of CVs in NGC6791 is consistent with 
a primordial population and does not require any dynamical 
formation or destruction processes. 

It is possible that there are more CVs in NGC6791 that 
we missed beca use they are to o faint in either X-rays or the 
optical. Byckling et al. (2010) studied a sample of twelve 
dwarf novae within ~200 pc with parallax distances. Only two 
(GW Lib and WZ Sge) have X-ray luminosities very close or 
below our detection limit, which is ~1 x 10 30 erg s (0.3-7 
keV) for the typically hard spectra of CVs. To our knowl- 
edge, there are no similar, published, studies of the X-ray lu- 
minosity function for magnetic CVs. Systems with moder- 
ate magnetic-field strengths (between 1-10 MG, the so-called 
intermediate polars) are thought to be generally brighter in 
X-rays than dwarf novae due to the higher mass-accretion 
rates; see e.g. Heinke et al. (2008). Polars (B > 10 MG) 
on the other hand can be faint when in the low state with 
L x .u « 10 29 - 10 30 e rg s" 1 . EU Cnc, the CV in M 67 ( [van 
den Berg et alj|2004[ ), and EFEri ( |Schwope et Hfffin) are 
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two such cases. While polars make up a relatively small frac- 



tion of known CVs (about 14% in the catalog by Ritter & Kolb 
(2003 )), their intrinsic contribution to less-biased CV samples 
is not well known. We estimate that a small (~15%) fraction 
of dwarf novae and an unknown fraction of polars could have 
gone undetected by our Chandra observation. 

The detection limit of the S03 catalog is V « 24, which 
corresponds to My w 10.5 at the distance of NGC 6791, 
taking into account extinction. We comp iled a sample o f 
27 CVs with parall a x distances fro m Patters on et al.|([2 008 ), 
Thorstensenj ( |2003) , |Harrison et al.| ( |20"0 4), and Thors tensen| 
et al.| ( |2008). Using the typic al magnitude of each system 
(MAGlln |Ritter & K6Tb] ( |2003) l), we find that the correspond- 
ing My distribution extends down to My * 12.2, with about 
25% of sources (including the X-ray faint systems GWLib 
and WZ Sge) below the S03 detection limit. 

While the above arguments sugg est that we should h ave de- 
tected most CVs in NGC 6791, [Gansicke et aij ( |2"00"9) found 
that the period distribution of Sloan Digitized Sky Survey 
CVs has a spike around the period minimum (80-86 min) 
corresponding to low-accretion rate systems. The incomplete- 
ness of our observations to such faint systems, which have an 
average magnitude of M g = 1 1.6 + 0.7, is significant. 

5.2. Active binaries in NGC 6791 

For ABs inside the half-mass radius of NGC 6791, the 
Chandra observation has a detection limit of about (0.7 - 1) x 
10 3() erg s _1 (0.3-7 keV) if their coronal temperatures are be- 
tween 1 and 2 keV. By comparison to the luminosity distribu- 
tion of ABs in M 67, this implies that many more ABs can be 
present in NGC 6791 than those listed in Table E] 

Fig. [7] shows the Lx versus orbital-period distribution for 
ABs in M67 and NG C 6791. The data for M 67 (open cir- 
cles) were taken from |van den Berg et aL] ( |2004| l, but, for 
easier comparison with the present analysis, the X-ray fluxes 
were recomputed such that they correspond to emission from 
a 2-keV plasma. The ABs in M 67 show a dominant trend of 
decreasing Lx with increasing period for orbital periods be- 
tween about 4-5 and 10 days. This can be explained by the 
activity-rotation relationship: stars in tidally-synchronized bi- 
naries with longer orbital periods rotate slower, and therefore 
are less active in X-rays, than stars in systems with shorter 
periods. For the shortest-period ABs, Lx goes down again, 
which is a result of saturated activity levels. The X-ray lumi- 
nosity of coronally-active stars do es no t exce ed 0.001 times 
the bolometric luminosity (e.g. |Giidel| ( |2004| l), which goes 
down for the shortest-period binaries. For contact binaries, 
the maximum X-ray luminosity is even lower, an effect known 
as super-saturation. In Fig. I7J stars that have evolved past 
the main-sequence turnoff or that are at the turnoff are plotted 
with larger symbols. In both clusters, these particular stars 
are only m oderately evol ved, with estimated radii no larger 
than ~3 R (Flower 1996), and all lie on or near the sub-giant 
branch. Except for two cases, the sub-giant ABs do not show 
a different trend than the main-sequence ABs (small sym- 
bols). The exceptions CX23 (V 100) in NGC 6791 and S 999 
in M 67 (upper right in Fig. [7) lie above the overall rotation- 
activity trend, and are good targets for follow-up study. 

The ABs in M67 show a clear separation in orbital pe- 
riods between detections and non-detec tions (not shown in 
Fig. [7] but see |van den Berg et al. ( 2004| l): practically all bi- 
naries with periods below the tidal-circularization cutoff pe- 
riod (12 d in M 67, |Latham et "aT1|1992| > are detected and, 
except for eccentric systems with pseudo-synchronous ro- 
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Figure 7. Lx versus orbital-period distribution of ABs inside the half-mass 
radii of M67 and NGC 6791. Small symbols are systems on the main- 
sequence, larger symbols mark systems at the turnoff or on the sub-giant 
branch. The eccentric binaries S 1242 and S 1272 in M67 are plotted at 
their pseudo -sync hronous rotation periods. The two possible periods for 
CX23 (Sect. [43) are connected with a horizontal line. Our sample of ABs 
in NGC 6791 (filled circles) is expected to be most incomplete for periods 
between a few days up to the circularization cutoff period, which agrees with 
the periods of the undetected periodic photometric variables in the cluster 
(filled triangles). The shaded area marks the upper limit for detection in 
NGC 679 1 . We hav e used the recent proper-motion membership study from 
|Yadav et al.| (2008) to re-assess membership for the M 67 Chandra sources 
and thus classified three more ABs with known periods as likel y cluster mem- 
bers; these a re CX 15, CX 58 and CX 61 (see Tables 2 and 3 in van den Berg 
|eTaLlp04) ). 

tation periods shorter than 12 d, all normal main-sequence 
and (sub)giant binaries with longer periods remain undetected 
down to Lx ~ (2 - 6) x 10 28 erg s _1 . This can be understood 
if the activity levels of the stars in these wider systems are 
closer to those of single stars in the cluster because their ro- 
tation is not (sufficiently) enhanced by tidal coupling. If the 
ABs in NGC 6791 from Table|2]show a similar trend as those 
in M67, it is clear from Fig.[7]that our Chandra observation 
of NGC 6791 can only see the brightest part of the AB pop- 
ulation, and that our AB sample is most incomplete for peri- 
ods ranging from a few days up to the circularization cutoff 
period in NGC 6791, which is at least 15 d. This is consis- 
tent with the period distribution of the detected ABs (filled 
circles) and non-detected periodic photometric variables with 
membership propability > 50% (filled triangles at the up- 
per limit of detection in Fig. 17). Variables for which the pe- 
riod is ambiguous are excluded from the figure, and so is the 
periodic variable V 78 which is mo re likely to be a pulsator 
rather than a binary (Kaluzny 2003). Not all ABs are neces- 
sarily photometrically variable all the time; a comprehensive 
radial-velocity survey or deeper X-ray data would be useful 
to uncover more of them. 

5.3. On the origin of sub-subgiants 

Sub-subgiants in star clusters are typically detected in X- 
rays, but the X-ray emission is explained by distinctly differ- 
ent types of binary interaction, making it a rather inho moge- 
nous group. Some resemble ABs ( see pVIathieu et al. (2003) 
for M67, and |Albrow e TaL] ( pOOl) and |Heinke et"ar| (2005) 
for 47 Tuc), others are CVs (AKO 9 and PC1-V1 1 in 47 Tuc 



Albro w~et al.| 2001) , and at least one but possibly two (in 
NGC 6397) are in binary MSPs where the X-rays originate 
in the shock region between the mass outflow from t he sec - 
ondary and the energetic pulsar wind (Bogd anov et al.|2010 l. 
The disturbance of hydrostatic equilibrium in the donor in re- 
sponse to mass transfer could explain the CMD location for 
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the CVs, but the sub-subgiants in M67 and NGC6791 show 
no hint of ongoin g mass transfer. Moreover, as has been noted 
by |Mathieu et a l. (2003), the eccentric orbit of S 1063 ar- 
gues against mass transfer in the past, as a (nearly) Roche- 
lobe filling star would have quickly circularized the orbit on 
a time scale of ~10 yr. The fact that we see at least five 
sub-subgiants in NGC6791 is an argument against invoking 
any short-lived perturbed state (resulting from either regular 
binary evolution or dynamical encounters with other cluster 
stars), and is more suggestive of a hitherto overlooked binary 
evolu t ionar y path. N-body simulations of M 67 by |Hurley| 
|et al.| ( |2005| l created one star below the subgiant branch that 
offers a possible clue. This simulated star formed when two 
stars in a binary merged after the donor ascended the sub-giant 
branch which resulted in a common-envelope phase. The 
merger product is a single giant with an undermassive core, 
which explains its location below the normal cluster sub-giant 
branch. If this scenario applies to the observed sub-subgiants, 
which all show signs of binarity, the systems should origi- 
nally be in triples. No primordial triples were included in 
the models by |Hurley et al.| ( [2005| l, but understanding the ori- 
gin of sub-subgiants would be a strong motivation for taking 
them into consideration in future cluster simulation s. As has 
been discussed i n the context of blue stragglers by Perets &| 
triples can indeed facilitate the formation 
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of very close inner binaries that eventually merge or enter a 
phase of mass transfer. 

5.4. X-ray source populations of old open clusters 

We compare the close-binary populations in three old open 
clusters— NGC 6791, M 67, NGC 6819— as probed by X-ray 
observations to look for trends with cluster properties, specif - 
ic ally mass or age. N GC 68 19 was recently studied in X-rays 
by |Gosnell et al.| ( [20T2] l using XMM-Newton data. As we dis- 
cuss in more detail below, compared to M 67 and NGC 679 1 , 
classification of the X-ray sources in the field of NGC 6819 is 
far less complete at this point, but we include the cluster as 
it extends the age range of the sample down to 2.4 Gyr. Pre- 
vious studies have shown that all three clusters have a high 
fraction of binaries: almost 40% in NGC 6819 (spectroscopic 
solar-type binaries; Math ieu|2008| l, -50% in M67, and >25- 
30% in NGC 6791 (photometric estimates from t he width o f 
the main sequence; Fan et al. 1996 and Bedin et al. 2008 ). 



In order to set a uniform sensitivity, we only consider X-ray 
sources brighter than Lx — 1 x 10 30 erg s . We do not in- 
clude clusters that have been observed with ROSAT only (see 
|Belloni| ( [l997[ ) for an overview) as these observations are gen- 
erally shallower, and the larger positional uncertainties result 
in less secure source identifications. 

Table PT] gives an overview of the cluster properties and 
the number of X-ray sources Nx, in total and separated by 
class, detected within the approximate half-mass radii. For 
NGC 6791 we refer to the list of sou rces (members and possi- 
ble members) in Table|2] Platais et al. (201 1) estimate a lower 
limit to the cluster mass of ~5000 M G ; taking into account 
detection incompleteness, stellar remnants, and binarity, the 
total mass is likely not higher than ~7000 M Q . For M 67 we 
adopt a tota l mass of -1100 M Q ( |Richer et aL]|1998| l. |Fan| 
|et al.| ( [T996| measured the half-mass radii for several types of 
stars (giants, single main-sequence stars, and main-sequence 
binaries); as a result of the effects of mass segregation, they 
find values ranging from ~7'2 to 12', depending on the typ- 
ical mass of the objects. For simplicity, we count the num- 
ber of Chandra sources inside the ACIS-I field (16'8 x 16'8), 



which is clo se to the area inside any of the r/, estimates by |Fan| 
[etaLl ( |T996l l. For both M 67 and NGC 6791 we do not count 
unidentified sources, or unclassified sources with candidate 
counterparts that lie away from the main sequence, as possi- 
ble members because the extra-galactic log N - log S curves 
sugg est they are most likely background AGN. 

|Gosnell et al.| ( |2012] l detect twelve sources inside r/, = 33 
down to L x = lxlIFergs- 1 (0.2-10 keV) in aXMM-Newton 
pointing of NGC 68 19. The estimated number of background 
sources suggests at most six to seven of these are cluster mem- 
bers. The five sources with candidate optical/UV counter- 
parts include potentially interesting systems like a candidate 
CV, a candidate qLMXB (the nature of both need to be con- 
firmed spectroscopically), a possible sub-subgiant (X 9), and 
an RS CVn binary that may have been formed in a dynam- 
ical encounter (X6). Given the depth of the optical catalog 
(V ~ 24) used to look for matches, it is reasonable to assume 
that an upper limit to the number of cluster ABs is the number 
of sources with a candidate counterpart on the (binary) main 
sequence or (sub-)giant branch that cannot be excluded as a 
non-member. There is at most one such system (X6). Due to 
the incomplete source classifications, the total cluster X-ray 
luminosity is highly uncertain, and we do not list it in TableB] 
If we simply take the sum of the luminosities of the twelve 
sources and divide it by two, we find log(2L3o/M) * 29.1 
(0.2-10 keV); as this estimate is already very approximate, 
we do not a ttempt to convert this value to the 0.3-7 keV band. 
|Kalirai et aL] ( |2001| l estimate the mass of NGC 68 19 to be 
~2600M o . 

The numbers of CVs {Nx,cv) an d sub-subgiants {Nx,s ) are 
smalQ only in NGC 6791 are there any spectroscopically 
confirmed CVs above our luminosity cutoff. A positive corre- 
lation of Nx.cv and Nx,s with cluster mass is allowed. Based 
on a compariso n wi th the field CV space density, we already 
argued in Sect. [57T] that the CVs in NGC 6791 could be pri- 
mordial, and the same is true for the (candidate) CVs in M 67 
and NGC 6819. This indeed implies that Nx,cv is proportional 
to cluster mass. 

The value of Nxab in Table HI gives the number of normal 
ABs, i.e. it does not include blue-straggler systems, which 
may have a non-standard evolutionary history, or the long- 
period AB-like sour ces (such as the yellow s tragglers S 1040 
and S 1072 in M 67; van den Berg et al. 2004} . Comparison of 
Nx,ab for the three clusters shows that a simple scaling with 
the present-day cluster mass is excluded. This is true for the 
ABs with evolved components (Nx,ab-g), which have wider 
orbits and may thus be easier involved in dynamical interac- 
tions, but also for ABs on the main sequence (Nx,ab-ms)- Al- 
though NGC 6791 is 4.5-6.4 times more massive than M67, 
it has 0.9-1.6 times the number of ABs. NGC 6819 has 2.4 
times the mass of M 67, but only has a fraction (~0. 13) of the 
number of ABs. These ratios show that Nx,ab normalized by 
cluster mass does not show a trend with age or mass either. 
M67 has the highest specific AB abundance (~0.014 per half 
the cluster mass), while this number is at least ~3 times lower 
in NGC 68 19 and NGC 6791, which are both more massive 
than M67, but younger and older, respectively 

It is unclear if the relative dearth of ABs in NGC 6819 and 



7 The CV EUCnc in M67 is fainter in X-rays, and strictly speaking it 
has not been demonstrated that it is a cluster member. The well-studied sub- 
subgiant S 1 1 13 in M 67 is not counted in TableHlas it lies 13'.7 away from the 
cluster center. Likewise, the sub-subgiants CX68 and CX 77 in NGC 6791 
lie at r > r/, . 
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NGC 6791 with respect to M67 can be explained by stellar 
or binary-evolution arguments, or if cluster dynamics plays a 
role. Possibly, M 67 is an outlier. Several studies hav e found 
that this cluster is dyn amically highly evolved (e.g. |Daven- 



port & Sandquist 2010). If, as a result of mass segregation, it 
lost a significant fraction of its mass due to the evaporation of 
preferentially low-mass stars, the current binary population 
would appear as representive of a much more massive clus- 
ter. To explain the discrepancy in the number of ABs, this 
effect should be more significant for M 67 than for NGC 6819 
and NGC 6791. Another difference is that the X-ray source 
population of M 67 contains binaries (blue and yellow strag- 
glers) with periods in excess of hundreds of days that look like 
coronal emitters (but note that these are not included in Nx,ab 
in Table [4]); such puzzling systems are not seen in the other 
two clusters. Clearly, a larger cluster sample spanning wider 
ranges in age and mass needs to be studied to make sense of 
what is the norm and what the exception, and understand the 
X-ray emission of old open clusters (to which ABs contribute 
an appreciable fractiorrji in general. As NGC 6791 is one of 
the oldest open clustersknown, we expand our comparison to 
globular clusters in the next section. 

There are two more old (> 2 Gyr) open clusters that have 
been observed in X-rays, but that we have not considered 
here. The young (1.7-2 Gyr) and nearby cluster NGC 752 
is spatially extended with members spread over a much larger 
region (~90' in diameter) than covered by the Chandra and 
XMM-Newton pointings. The census of X-ray sources could 
thus be very i ncomplete. We note that the list of detected 
members from Giardino et al. ( 2008| ) includes just one source 



with Lx > 1 X 10 3 " erg s , which is matched to a likely 



main-sequence AB. jGondoin ( 2005| > analyzed XMM-Newton 
data of NGC 188 with a sensitivity of L x = 9 x 10 29 erg s _1 
(0.5-2 keV), about six times deeper than the ROSAT study by 



Belloni et al.| ( |1998) . However, the chosen aimpoint is offset 
from the cluster center by almost 6' . Only new observations 
can tell whether this explains why no sources are detected in 
1/4-1/3 of the contiguous area inside the half-mass radius that 
is farthest from the observation aimpoint. In any case, this 
unfortunate circumstance makes NGC 188 less suitable to be 
included in our comparison^] 

5.5. Comparison with globular clusters 

While almost 80 globular clusters have been studied with 
Chandra, only a handful have been observed deep enough 
to uncover a significant fraction of the ABs. Here we focus 
on two globular clusters for which the Chandra sensitivity 
reaches well below the luminosity cutoff adopted in the previ- 
ous section. The nearby core-collapsed cluster NGC 6397 was 
observed down to a lim iting luminosity of ~9 x 10 28 erg s _1 
(Bogd anov et al.|2010|i. We a d opt a cluster mass o f 2.5 x 10 5 
Mq ( [Pryor & Meylan|[l993) . |Cohn et al.| (pOKj) clas sified 
most of the 79 Ch andra sources inside r/, - 2'. 33 from Bog- 
|danov et aL] ( |2010| using H5Tdata, and found numerous CVs 

8 -15% in NGC 6791, -45% inM67 

9 Moreover, some sources seem to be missing from the list of six detected 
cluster members in Gondoin 1 2005 1. ROSAT source X29 in [Belloni et al.l 
(1998) , which is the brightest cluster X-ray source, is not included in the 
source list despite being clearly visible in the EPIC-MOS images; in the 
XMM-Newton Serendipitous Source Catalog (Watson et al. 20091 it has a 
maxim um source de tection likelihood of 279. ROSA1 source X 21 is identi- 
fied b y|Belloni et al.| t 19981 with a likely short-period AB and cluster member. 
The off set between this optical source and the corresponding XMM-Newton 
source S 17 is 3.2", and it is unclear why it was not identifie d as a potentia l 
counterpart given the match radius of 8" that is adopted by Gondoin| (2005) . 



and ABs, and two sub-subgiants (at least one is in a binary 
with a MSP). About 20% of the sources are brighter than our 
cutoff value; these are mostly CVs or neutron-star binaries. 
While there are plenty of X-ray faint ABs, there are just a few 
at most above 1 x 10 30 erg s" 1 ; the range listed in Table |4] 
is zero to two, but those two systems have very red colors 
and are unlik ely cluster members. 47 Tuc is a massive cluster 
(1.3xl0 6 M Q ; |Pryor & Meylan|l993| . We estimate that within 
r/, = 2!79, about 28 of the 232 sources bright enough to be in- 
cluded in our selection are background AGN, leaving ~200 
sources associated with t he cluster. As of yet, o nly part of 
these have been classified (Edmonds et al. 2003a b). We make 
a very conservative estimate of the number of CVs and ABs. 
A lower limit is derived by cou nting the n umber of (candidate) 
CVs and ABs in Table 2 of Heinke et al.| ( |2005| l, while an up- 
per limit (given in square brackets in Table [4]i is estimated by 
adding all unidentified and unclassified sources to either the 
CVs or ABs and subtracting the estimated contribution from 
AGN. Sub-subgiants are more conspicious, so their number is 
less uncertain. To compute the cluster X-ray luminosity, we 
sum the luminosities for the Nx individual sources, and sub- 
tract an approximate contribution for the background AGN 
based on the mean luminosity of the unidenti fied sources. For 
both c luste rs we converted the fl uxes listed in Bogdano v et al.| 
(2010) and jHeinke et al. (2005) to the 0.3-7 keV band using 
PIMMS, adopting a 2-keV plasma model. 

Table H] readily shows that the total number of ABs above 
Lx = 1 X 10 30 erg s _1 , normalized by cluster mass, is 
much lower in the two globular clusters than in M67 and 
NGC 6791. For 47 Tuc the difference is a factor of 65-230 
compared to M 67, whereas NGC 6397 likely does not have 
any ABs at all above this luminosity. Lowering the X-ray 
luminosity cutoff value would increase the numbers of ABs 
especially for NGC 6397, but this is not enough to make up 
for the lower specific X-ray luminosity. The poor statistics in 
the open clusters left aside, CVs appear underrepresented in 
globular clusters as well, by about an order of magnitude. A 
likely explanation is that dynamical encounters have a net ef- 
fect of destroying ABs — both those with evolved components 
and wider orbits and main-sequence systems — and CVs. At 
the same time, the c orrelation with encounter rate found by 
|Pooley & Hut| ( |2006|l implies that at least CVs are formed 
as well; perhaps this is why the specific frequencies in open 
and globular clusters are less disparate for CVs than for ABs. 
AB-like sub-subgiants are underabundant in globular clusters 
as well, but on the other hand sub-subgiants in (candidate) bi- 
nary MSPs are not found in open clusters (note that in Table|4] 
the different types of sub-subgiants are grouped together be- 
cause some have not been studied well enough to make the 
distinction). 

The relative lack of ABs in NGC 6397 with respect to other 
globular clusters was already noted by Bassa et al. (2004]) and 
Grindlay (2006), who ascribed it to its post-collapsed state. 
We now s how that 47 Tuc is deficient in ABs, too. M4, 
studied by Bassa et al. (2004) and not included in Table [4] 
shows the same effect. Even if all fifteen sources or so that 
are brighter than 1 x 10 30 erg s _I were ABs (which is not 
the case as there are also CVs, background and foreground 
sources, and an MSP among them), it is not enough to reach 
the number of ABs expected based on scaling by mass from 
open clusters. Milone et al. ( 2012) ) find a lower ov erall bi 



nary fractio n in globulars compared to open clusters (Sollima 
et al.|2010 i. If this is the result of binary destruction, the X- 
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Table 4 

X-ray sources in globular and old open clusters (r < r/,, L\ ^ 1 X 10 30 erg s ) 



cluster 


age 


M 


N x 


Nx.cv 


N x ,s 


Nx,AB 


NxAB-MS 


NxAB-G 


log(2 L 30 /M) 




(Gyr) 


(M e ) 
















NGC6819 


2-2.4 


2600 


6-7 


J? 


J? 


1? 





1? 




M67 


4 


1100 


12 





1 


7-8 


5-6 


2 


28.9 


NGC6791 


8 


5000-7000 


15-19 


3^ 


3 


7-11 


5-9 


2 


28.6-28.8 


47Tuc 


11.2 


1.3 X 10 6 


-200 


30[-119] 


10 


42[-131] 


34 


8 


28.0 


NGC 6397 


13.9 


2.5 x 10 5 


15-18 


11 


2 


0-2 


0-2 





27.7 



Note. — We list the number of sources identified with cluster members (Nx), CVs (Nx,cv), sub-subgiants (Nx,s), and ABs (A'x.ab)- We sub-divide the ABs 
into those on the main sequence (Nxab-ms) and evolved systems (Nxab-g)- The numbers for NGC 6819, printed in italics, are especially uncertain due to the 
limited optical follow up. Log (2 L30/M) is, in log units, the ratio of the total X-ray luminosity of the Nx sources with Lx <; 1 x 10 30 erg s - ', to the cluster mass 
(M) divided by two to account for the selection of sources inside the half-mass radius. The X-ray luminosity refers t o the 0.3-7 keV band and corresponds to 
emission from a 2 keV plasma. The horizontal li ne sep arat es the open c lusters from the globular clusters. References: |Gosnell et al.| (2012 1 for NGC 6819, [van] 
|den Berg et aT](2004) for M 67, |Bogdanov et al.|j2010| and |Cohn et aTlpOlO) for NGC 6397, |Heinke et al.|j2003) for 47 Tuc. 



ray emitting, and thus hardest, binaries are affected as well. 
We cannot exclude that other factors are important. There is a 
trend of decreasing AB frequency with age for all clusters in 
Table [4|excep t NGC 6819. 

In |Verbunt| d2000) it was first pointed out that the total 
ROSAT luminosity per unit mass of most globular clusters is 
lower than that of M 67. While M 67 may be exceptional in 
some ways (see Sect. 5.4 1, the last column of Table [4] shows 
that the specific X-ray luminosity of globular clusters may 
be lower than for open clusters in general. CVs, ABs, and 
sub-subgiants are all responsible for this. The X-ray emission 
from qLMXBs and MSPs in globular clusters, of which there 
are no confirmed cases in old open clusters, cannot make up 
for the lack of other types of faint sources. A more extensive 
investigation, based on a larger sample of globular clusters 
and old populations, is done by Heinke et al. (in preparation) 
to investigate the effects of age, central density, and metallic- 
ity on differences in the X-ray emissivity. 

6. CONCLUSIONS 

We have studied NGC 6791, which is one of the oldest open 
clusters known, in X-rays for the first time. Of the twenty 
Chandra sources that are identified with proper-motion clus- 
ter members, nineteen are binaries. We find a mix of CVs, 
ABs, and sub-subgiant binaries. With optical follow-up spec- 
troscopy we confirm the classification of one CV in the clus- 
ter, and identify one new candidate CV. The origin of these 
and the other two cluster CVs is likely primordial. A compar- 
ison with the X-ray sources in the younger cluster NGC 6819 
suggests that in order to find numerous (main-sequence) ABs 
brighter than Lx — 1 x 10 30 erg s _1 in open clusters, one 
must look at those older than ~2.5 Gyr. The mass-specific 
frequency of main-sequence and sub-giant ABs in NGC 6791 
(8 Gyr) is 3-7 times lower than in M 67 (4 Gyr); it is not clear 
why this is the case. Comparison with a small number of glob- 
ular clusters shows that all three source classes that are mainly 
responsible for the X-ray emission from old open clusters 
are under-represented in the globulars. This accounts for the 
lower total X-ray luminosity per unit mass of globular clus- 
ters, and indicates that the net effect of dynamical encounters 
may be destruction of-even the hardest-binaries. The cutoff 
luminosity for comparing source populations among clusters 
wassetatLx = lxlO 30 erg s -1 (0.3-7 keV). A lower cutoff lu- 
minosity might be needed to account for effects of variations 
of the AB X-ray lum inosity function with age or metallicity. 
While Verbunt ( 2000[ > note that population-II ABs may have 
lower X-ray luminosities, it remains unclear in this context 
why NGC 6791, which has a super-solar metallicity, should 



also have relatively fewer ABs above our limiting L x . We de- 
fer a more thorough investigation to a future study. Deeper 
data are needed for NGC 6791 to match the already deep ex- 
isting data for M67, 47 Tuc, and NGC 6397. The X-ray de- 
rived AB populations may then be compared to those derived 
from other wavelengths, e.g. optical variability studies. Also 
in this respect, studies of NGC 679 1 may benefit from being 
in the field of view of the Kepler satellite. 
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facility of the Smithsonian Institution and the University of 
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